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Abstract 

We derive two classes of brane-world solutions arising in the presence of a bulk scalar 
field. For static field configurations, we adopt a time-dependent, factorizable metric ansatz 
that allows for radion stabilization. The solutions are characterized by a non-trivial warping 
along the extra dimension, even in the case of a vanishing bulk cosmological constant, and 
lead to a variety of infiationary, time- dependent solutions of the 3D scale factor on the 
brane. We also derive the constraints necessary for the stability of these solutions under 
time-dependent perturbations of the radion field, and we demonstrate the existence of 
phenomenologically interesting, stable solutions with a positive cosmological constant on 
the brane. 



1 Introduction 



Over the last few years, there has been considerable interest in models in which our universe 
is a 3-brane (a hyper-surface) embedded in a higher dimensional bulk. Much of the interest 
in extra-dimensional field theory is due to the hope for a solution to the hierarchy problem 
[1, 2, 3]. In these models, the extra dimensions are hidden from us, not necessarily by their 
smallness but by our confinement to a four- dimensional slice of the bulk spacetime [1]. In 
contrast to Kaluza-Klein scenarios, standard model interactions are confined to a brane 
whereas gravity propagates through the bulk perpendicular to the brane. The hierarchy 
problem can be resolved by either postulating large extra dimensions (in which case the 
TeV scale is the fundamental scale of gravity and the Planck scale is derived in terms of 
the fundamental scale and the volume of the extra-dimensional space) [1] or when the 4D 
metric scales exponentially throughout the bulk (the so called "warp" factor) [2]. 

While static brane-world models have served as a useful tool for testing ideas in higher 
dimensional spacetimes, their direct applicability to cosmology is limited. More realistic 
cosmological models may be derived by allowing a non-vanishing four- dimensional cosmo- 
logical constant, or by introducing time-dependent energy-densities on the branes. Various 
cosmological aspects of such models have been investigated in the literature [4]-[12]. One 
of the serious problems in brane models is the resulting unconventional set of Priedmann 
equations [5, 6, 7]. The Hubble parameter, H, on the brane is often found to scale as H p, 
rather than the standard four dimensional dependence, H ~ ^/p. It has been shown how- 
ever, that this problem can be solved upon the proper stabilization of the extra dimension 
[13]-[15], that removes any unnecessary constraints between the brane energy-densities. 

Both static and time-dependent generalizations of the original Randall- Sundrum solu- 
tions [2] have been also constructed by introducing a bulk scalar field [16, 17]. As a matter 
of fact, the task of the stabilization of the extra dimension was first accomplished by in- 
troducing a bulk scalar field, which had different non-vanishing vacuum expectation values 
on each of the two branes [18]. The same topic was further elaborated in [17, 19]. Here, 
we will derive two classes of brane-world solutions which include a static bulk scalar field 
and provide infiationary solutions in the 4D slices. The two classes correspond to either 
vanishing or non-vanishing bulk potential for the scalar field, and are characterized by a 
non-trivial warping of the metric along the extra dimension even in the case of zero bulk cos- 
mological constant. Due to the appearance of a bulk curvature singularity, we are forced to 
consider two-brane-system configurations which, however, have a fixed inter-brane distance 
and lead to conventional FRW equations on the branes without any additional fine-tuning. 
By using a method developed recently in Ref. [20], we study the stability of those solutions 
under time-dependent perturbations of the radion field, and demonstrate that we can easily 
find parameter regimes where these solutions are stable. Even more important is the fact 
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that some of these stable solutions have a positive cosmological constant on the brane - 
previously all known solutions of this type were unstable [21]. 

We organize this paper as follows. In section 2, we present the equations of motion of 
our theory and show how a factorizable (in time and the extra space coordinate y) scale 
factor can be obtained in the presence of a bulk scalar field. In section 3, we present 
two classes of inflationary brane solutions involving a static bulk field with vanishing or 
non-vanishing, respectively, bulk potential and we demonstrate that these solutions indeed 
lead to conventional FRW equations. We derive the conditions for the stabilization of 
these solutions in section 4 and investigate the parameter regimes that correspond to stable 
configurations. Finally, in section 5, we summarize our results. 



2 Equations of Motion for Gravity and a Bulk Scalar 
Field 

We start from the cosmological principle of isotropy and homogeneity in the three space- 
like dimensions of the brane. The presence of the brane breaks the isotropy along the fifth 
dimension and this is refiected in the explicit y-dependence of the metric tensor. Based on 
these facts, we make the following ansatz: 

ds'^ = Qmn dx^dx^ = -ri^it, y) df + a^{t, y) 7^^- dx'dx^ + 6^(i, y) dy"^ , (2.1) 

where M,N = 0, 1,2,3,5, 'jij is the usual Robertson- Walker 3-space metric tensor, and t, 
(i = 1, 2, 3), and y are the time- and space-like coordinates along the brane and the extra 
dimension, respectively. 

In addition, we consider a bulk scalar field (f){t, y), which depends only on time and the 
extra coordinate^. The action of this five-dimensional, gravitational theory is given by: 

r ( - 1 

S=- J d'xdy^l—^R + AB + -dM<l>d"4> + VB(<l>) 

+ 5:[A. + V^W]^|. (2.2) 

i ) 

where M5 is the fundamental, five-dimensional Planck mass, R denotes the five-dimensional 
scalar curvature, and Vb and Vi stand for the bulk and brane potentials, respectively, of 
the scalar field. Finally, and Aj are the vacuum energies of the bulk and the branes. 
From Eq. (2.2), one can derive the scalar field equation of motion: 

^ Again, we use the standard assumption that the 3-dimensional space is homogeneous and isotropic. 
Thus, the field is independent of the 3-dimensional spatial variables. 
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where dots and primes denote derivatives with respect to t and respectively. 

The matter content of the five-dimensional space-time is described by the energy- 
momentum tensor of the bulk scalar field and the bulk cosmological constant, which may 
be written as: 



(2.4) 



where /i,!/ — 0, 1, 2, 3, and 

Tmn{(I)) = (?M0 (?JV0 - Qmn 



(2.5) 



We next consider the five dimensional set of Einstein equations. They can be obtained 
by the variation of the action (2.2) with respect to the metric, and, for the metric ansatz 
(2.1), they take the form ^ 
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SttGs = Sir/M^ is the five- dimensional Newton's constant, and k 
denotes the spatial curvature of the four- dimensional spacetime along the brane. 

While the metric is continuous across the extra dimension, its derivatives with respect 
to y can be discontinuous because of the inhomogeneity of the matter distribution in the 
fifth dimension, notably the branes. Therefore, a delta function appears in the Einstein 
tensor and this must be matched with a delta function in the energy-momentum tensor 



^Our convention for the Riemann curvature tensor is 
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[22]. Here we use a similar notation for the jump of the scale factors^ as in references [7, 13] 
and obtain 
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where the subscript i denotes evaluation of all quantities at y = yi. The corresponding 
jump condition for the scalar field (p is obtained from the field equation (2.3) 
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(2.12) 



y=yi 



While the jump conditions for the scale factors (2.10) and (2.11) are non-trivial given any 
inhomogeneous source on the brane, the jump condition for the bulk field (2.12) depends 
only on the interaction between the bulk field and the brane. In particular, if the bulk field 
is sitting at an extremum on the brane, then the jump of the scalar field (2.12) is zero. 

The facility to find an exact solution to the field equations is often aided by our ability 
to factorize the scale factor on the brane. Consider the following. 
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If the scale factor can be factorized, then the right-hand-side of Eq. (2.13) must vanish. 
From Eq. (2.13), we note the following: if the field is sitting at a local minimum on the brane 
or the field is static on the brane, our ability to factorize the scale factor into independent 
t and y dependencies is tied to the stability of the extra dimension on the brane {pi = 0). 
Off the brane, a corresponding relation can be found by examining the Gqs component of 
the Einstein equation (2.8), which can be rewritten as 



n' a' \ a, a' b d j a' 
n a I a ah dt\ a , 



(2.14) 



Without loss of generality we can write n{t, y) — n{y). Then, Eqs. (2.10)-(2.11) lead to the 
following factorizable form of the 3D scale factor ^ 



a{t,y) = a{t)n{y) 



(2.15) 



^Here a" = a'ji+ [a']i 5{y — yi) where a'^ is the non-distributional part of the second derivative of a, and 
[a']i is the jump in the first derivative across y = yi, defined by [a']i = a'{t,y = yi + e) — a'{t,y = yi — e). 

"^A static four-dimensional universe can be obtained by n{t,y) = a{t,y) = n{y). This was the ansatz 
used in Ref. [2] to solve the hierarchy problem and obtain conventional Newtonian gravity. 



4 



Inserting the above into Eq. (2.14) shows us that 



no 6 



(2.16) 



Thus, if the bulk field is either time-independent or y-independent, our condition for the 
factorization of the scale factor is intimately connected to the stability of the fifth coordi- 
nate. In what follows, we are going to assume that = 6 = 0, which brings our metric 
ansatz to the form 



ds^^n^{y) -df + a^{t)-fijdx'dx^ +b\y)dy'^ 



(2.17) 



3 Exact Time-Dependent Solutions 

We now proceed to derive exact solutions of the coupled system of gravitational and scalar 
field equations that have a non-static, four-dimensional line-element. We are going to 
present two classes of solutions: the first class corresponds to a vanishing potential, Vb{(I)), 
for the bulk scalar field, while the second one arises in the presence of a non-trivial (expo- 
nential, in a certain limit) bulk potential. 



3.1 Solutions with a vanishing bulk potential: Vb = 

In our first example, the assumption of Vb(0) = leads to a simphfication of the field equa- 
tions in the bulk. The derivation of the solution is further facilitated by a transformation 
of the y-coordinate that allows us to write the metric ansatz (2.17) in terms of "conformal" 
coordinates, i.e. b{y) — n{y). Then, Einstein's equations may be rewritten as (henceforth, 
a denotes the 3D scale factor appearing in the metric ansatz (2.17)): 

3n" fa" ^^'^ 



n 



3n" /2a k\ o /^'^ \ /o 

- + 3 + 7^ H-'^U^+^^fi ' (3-2) 



n Va 'a^ a'^ J ^V2 
6n'2 d k- --^'2 



By combining Eqs. (3.1) and (3.2) we obtain 



-3K + 7 + 72 =^^5 V-^'As . (3.3) 
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from which we can find the following time-dependent, infiationary solutions for the scale 
factor ^ 

jj sinh[H{t - to)] when k ^ -1 (3.5) 

^ cosh[H{t - to)] when k ^ +1 

where H and to are constants. For A; = 0, we may obtain Minkowski, de Sitter or Anti 
de Sitter solutions on the brane when has either a zero, positive or negative value, 
respectively. For k = —1, we may have solutions for either positive or negative values of 
(in the latter case, the sinh-like solution above is replaced by sm[H{t — to)]/H, with 
= —H^ > 0) while for /c = +1 only solutions with > are allowed. 

Turning to the scalar field equation of motion (2.3), we see that it now becomes 




:0. (3.6) 
The above equation can be integrated once to give 

where c is an integration constant. The relation between (f)' and n, derived above, can be 
used in Einstein's equations to derive the form of the warp factor. In addition, we note 
that all solutions for the scale factor a{t) appearing in Eq. (3.5) satisfy the relation 

-= - +--H\ 3.8 
a \a J 

Inserting both (3.7) and (3.8) into the simplified set of field equations (3.1)-(3.3), we obtain 

,2 





+ n'ABj. (3.10) 

The above equations can be solved if we set = 0. In that case, we find the following 
solution: 

n^(y) = ^^^Hmy\) ^ H' = ^c'smh\3H\y^\), (3.11) 

sinh(3/J|yi|) ' 12 ^ ' ^ ' 

for a four-dimensional de Sitter spacetime {H^ > 0). The corresponding solutions for Min- 
kowski and Anti de Sitter 4D spacetime are straightforward to derive by taking the limits 

^In order to address the graceful exit problem, we would be required to consider a non-static bulk field 
which is beyond the scope of this paper. 
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H ^ and H — > iH, respectively. In the former case, the warp factor is hnear in y, 
while, in the latter, the solution is given in terms of a sin-like function. For simplicity, we 
will, for the remainder of this section, concentrate on the case with > 0. It is worth 
noting that, although we have assumed a vanishing bulk cosmological constant, non-trivial 
warping arises in all cases where the expansion rate, H, on the brane is non-zero. As can 
be seen from the second of Eqs. (3.11), the expansion rate is closely related to the kinetic 
term of the scalar field in the bulk which, in a way, replaces the bulk cosmological constant. 

The Ricci scalar of the five-dimensional spacetime described by the metric ansatz (2.17) 
has the form : 

^ 1 /6d^ 6a 6A;\ /4n'2 8n"\ ,^ 



The second equality can be obtained by using Eqs. (3.2)-(3.3), together with (3.4), with 
Ab = 0, in Eq. (3.12). As we can see, there is a true singularity in the bulk, at |/ = 0. If 
we choose to place the first brane at y — yi < (in which case, we have the normalization 
^(2/1) — 1)) then, in order to remove this singularity, a second brane must be introduced at 
a point y — y2, with yi < ^2 < 0, so that the singular point is never encountered. 

We next examine the jump conditions imposed on the warp factor and the scalar field 
at the positions of the two branes at yi and y2, respectively. Starting from Eqs. (2.11) and 
(2.12), using the 'conformal gauge' b{y) — n{y), and inserting the solutions found above, 
these conditions can be rewritten as 



y, + A,^^coth(3i/|y,|), 2c- ^^^^^^ 



kI d(j) 



y=yi 



(3.13) 



(3.14) 



y=y2 

where we have used the normalization condition ni = 1. As in the case of the two-brane 
static solution with a bulk scalar field found in Ref. [17], the form of the interaction terms, 
Vi, completely determines the ratio of the warp factors evaluated on the branes 

4 

/nix4 _ / Sinh(3gyi) y ^ idcl>V2)y=y2 /g^.^^ 

^n2^ \smh{3Hy2)) {d^V,)y=y,' ^' ^ 

Here too, the derivatives of the interactions are required to have opposite signs, just like the 
total self-energies of the two branes. Last, but not least, we may also observe that the inter- 
brane distance can be derived from the above relation which is accompanied, through the 
warp factor jump conditions, by the fine-tuning of only one of the two brane self-energies. 
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Let us briefly comment on the relation between the expansion rate H and the brane 
potentials VJ. To do so, we must first define the four-dimensional Newton's constant. Recall- 
ing the gravitational part of our original action (2.2) and noting that the first combination 
inside brackets in Eq. (3.12) stands for the 4D scalar curvature, we may write 

^Jd'xdy^R = ^Jd'xdy^,n'{^n^^^ + ---) 

= ^Jd'xV^^n^'^ + ■■■), (3.16) 
where ^—g/^ — a^{t). Thus, we have 



1 1 f 2 /■J'2 2 TV^ 
— = — fdyn^ = ^ dyn^^ ^ -j- coth{3H\yi\) 

J J y\ ACg Oil 



y2 



3/1 *=1 



where we have used the jump conditions (3.13) and (3.14) to obtain the last equality. Note 
that although it appears that is independent of ^5, {Vi + A^) and H are interdependent 
through K5. 

Finally, to check the Hubble equation on the brane, we must compute the effective 
cosmological constant, Ag//, defined as 

Aeff^fdyn'^-^ + Lj, (3.18) 

where M^'' is given by the second two terms of Eq. (3.12) and £ contains all of the non- 
gravitational pieces of (2.2). The result of the integration (taking care to eliminate the 
boundary terms associated with n") is 

Ae// = Enf^^, (3.19) 

i=i 

which when combined with Eq. (3.17), leads to the standard form of the Priedmann equa- 
tion, 

^eff 



H'-^Aeff (3.20) 



without the need of any fine-tunings. 

We will return to this solution in section 4, to test its stability under small time- 
dependent perturbations. 

3.2 Solutions with a non-trivial bulk potential: Vb 7^ 

In this subsection, we allow for a non-zero bulk potential, Vb ^ 0, a fact which will modify 
the equations in the bulk and subsequently their solution. Although we retain the form of 
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the metric (2.17), we choose to work with non-conformal coordinates; therefore, we redefine 
the y-coordinate and fix the scale factor 6 to a constant value, i.e. 6 = 1. Under the above 
assumptions, Einstein's equations become 

-3n^ (- + K]+^(-, + \]=4nHC + yB + AB), (3.21) 



\ n V? J \a? a? / % ^ 

( V + 7^) - (t + ^ + ^) - n\^^V.^Ks], (3.22) 

77/2 ,^2 a /C\ 2 2 



— - 3 - + - + - = - Vb - As . (3.23) 



By combining Eqs. (3.21) and (3.22), we can once again derive Eq. (3.4) indicating 
that our previous solutions for the scale factor a{t), appearing in Eq. (3.5), are still valid. 
What remains to be found is the new form of the warp factor and the scalar field. Plugging 
in the solutions for a(t), the above gravitational equations become 



2„- 2„^ |f_^^V,,(,) + ^,V (3.25) 



n 



3 V 2 

For H — Vb — 0, the above equations lead to the static solutions with a bulk scalar field 
found in Ref. [17]. In general, Eqs. (3.24)-(3.25) are difficult to solve, nevertheless, we can 
find a sub-class of solutions which satisfy the following: 

^ + Ivb{cI>)^E, (3.26) 

where £^ is a constant. The above relation can be used to provide the solution for both the 
scalar field and the warp factor. In the bulk, the scalar field equation of motion has the 
form 

By differentiating Eq. (3.26) with respect to y, we get 

^^0' = O, (3.28) 

which, when combined with Eq. (3.27), yields the simple equation 

n' 

0" + -0' = O, (3.29) 
n 
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whose solution is 

<t>'{y) = ^ , (3.30) 
niy) 

where c is again an integration constant. 

Going back to Eqs. (3.24)-(3.25), we can now combine them to get a single differential 
equation for niy) 

- = -^(3£; + Ab), (3.31) 
n D 

leading to the following solution for the warp factor 

where x^vl = 1 , 
'^^^^^ where sm^{ujyi) = oo^ , (3.32) 



n{y) 



sm(ui\yi\) 

smh{w\y\) 
smh{u>\yi\) 

In the above, we have imposed the normalization condition rii — n{yi) = 1, and defined 



^'"^(^1^1) where sinh' {uy^) = 



u;'^-^{3E + Ab) (3.33) 

D 



and 



X' = H'+'^. (3.34) 

The three solutions presented in Eq. (3.32) correspond to the combination (3£' + A^) 
being zero, positive or negative, respectively. All three solutions are characterized by a true 
singularity at y = 0. If we place, as in the previous section, a brane at y = yi < 0, then, 
a second branc should be introduced at y = ^2, with < 7/2 < 0, if we want the singular 
point shielded. 

The bulk potential of the scalar field is defined through Eq. (3.26). The exact form can 
be easily derived by using the expression of (f)', in terms of n{y), according to Eq. (3.30). 
The potential is everywhere well defined apart from the regime close to the singularity. 
Near y = 0, all three solutions (3.32) lead to the following expression for the bulk potential 

VB{<f>) = 3£; - I c'yl exp(^) . (3.35) 

^ ^\yi\ 

The scalar field, near the singularity, behaves as: (f) ~ — c|?/i| In \y\, and therefore diverges, 
as ?/ — 0, causing the bulk potential to diverge as well. The sign of c^, appearing in front of 
the exponential, determines whether the potential diverges towards plus or minus infinity. 
For > 0, 0'^ is also positive and the kinetic term of the scalar field has the correct sign; 
however, the potential diverges towards minus infinity, being unbounded from below. In 
this case, the introduction of a second brane to shield the singular point in the bulk is 
imperative. If, on the other hand, we consider the < case, then, we end up with a 
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'tachyonic' kinetic term for the scalar field; nevertheless, an infinitely-high potential barrier 
is rising in front of the singularity in this case, thus shielding the singular point and allowing, 
in principle, single brane configurations. Nevertheless, in what follows, we will introduce a 
second brane to shield the singular point even when the latter case is considered. 

The jump conditions, for the solution presented above, follow from Eqs. (2.11) and 
(2.12) upon substituting the expressions for the warp factor (3.32) and the scalar field 
(3.30) and setting b = 1 (here, we show the result only for the case 3E + Ab < 0). They 
take the form 



Fi + Ai = ^ coth(a;|yi|), 2c- '^^^ 



«5 



(3.36) 

y=yi 



V2 + A2^ coth(a;|2/2|) , 2c = -n2 



4 V".-./> 



(3.37) 



y=y2 

and may lead to the fixing of the location of the two branes in terms of the remaining 
fundamental parameters of the theory. 

Let us finally check the form of the Friedmann equation on the brane in this case. As 
in the previous subsection, we must first define the 4D gravitational constant. In terms of 
non-conformal coordinates, Eq. (3.16) still holds, with the factor n^, appearing in the first 
line, replaced by n'^. Then, we may write 

;^ = ^/^^" =^^(^^'^^)=.isinh^(c.,,)^S"'"^' ^ ^ ^ 

where we have used again the solution for cj^ > 0, and defined X(i/i, 1/2) as the integral of 
over the internal compact space. Turning to the form of the effective cosmological constant 
on the brane and using the fact that now Riy) = -(l2n'Vn2 + 8n"/n), we find 

Kff^^Ayi:y2). (3.39) 

In the above, we have also used the jump conditions (3.36)-(3.37). Eliminating the integral 
2^(1/1,1/2) from the above equation by making use of the definition of the 4D Newton's 
constant (3.38), we recover once again the conventional 4D Friedmann equation (3.20). 



4 Stability analysis 

In this section, we perform a stability analysis of the solutions found in the previous section 
under small, time-dependent perturbations by using the stabilization constraints that were 
recently found in [20]. Here, we will closely follow the method and notations used in that 
work deviating only when the needs of the particular solutions presented in this article 
demand it. 
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4.1 Stability behaviour of the solutions with Vb = 



We start our analysis by perturbing the metric ansatz (2.17), written in terms of conformal 
coordinates, in the following way 

ds'^ = n\t, y) [-dt^ + a^{t) dx'dx^ + }?{t) dy^\ . (4.1) 

As in [20], we are assume that the time-dependence of the scale factor h along the extra 
dimension induces a time-dependence to the warp factor n. Although the method that we 
follow is identical to that in [20], we need to repeat part of the calculation: a different 
metric ansatz is used here, a fact that changes the final expressions of the stabilization 
constraints. 

Starting from the gravitational part of the action, we need to express the five-dimensional 
scalar curvature in terms of the metric function appearing in (4.1). It takes the form 

R = ^7^(") - \d^D^% + R^y^ + Lun , (4.2) 


where 

In Eq. (4.2), Ckin contains time-derivatives of n and b that will eventually give rise to 
the kinetic term of the canonically normalized radion field. Here, we are interested in the 
stability behaviour of the solutions which can be found from the expression of the effective 
potential, therefore, we highlight only the points of our analysis that lead to the form of 
this quantity. 

Now, the gravitational effective action can be written as : 

Sg = -J d''xy/^^\^-A{b)'^ - dyn%{Ckin + R^"^)] (4.4) 

where ^J—g/^ — a^, and the conformal factor ^4(6) is : 

A{h)^%idyhn\t,y). (4.5) 

A conformal transformation of the four-dimensional metric, ((74)^^ = {cji^^yj A{\)\ removes 
the coupling between h and VS^^ and brings the total effective action to the form 

S^-\ ^i'x^/=|I{-^ + Ain + V;//}, (4.6) 
where the effective potential K// is given by: 



12 



The integral over y is performed over a compact dimension, with y taking values in the 
'circle' consisting of the symmetric intervals (^1,^2) and (y2,yi) while passing across the 
branes. 

The extremization and stabilization constraints follow from Eq. (4.7) by taking deriva- 
tives with respect to b. Wc may use the fact that b is always multiplied by y in the action 
to define a new coordinate ^ = by. Then, the effective potential may be rewritten as 



A\byi,by2)Veff 



/•by 2 

/ 

Jbyi 



1 - 
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k'2 



^2^ + ^)+i:Aby^){v. + A.), (4.8) 



where now primes denote derivatives with respect to ^. The extremization constraint follows 
by taking the first derivative of (4.8) with respect to b. Then, we obtain: 



eff 
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1=1 V 
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In order to simplify the above expression, we may now use the jump conditions in 
conformal coordinates, that relate and 0^ to V^'s and which are vahd for the static 
solutions. They can be written as 



rii 6 



, ni dVi{(l)) 



(4.10) 



y=yi 



for i = 1,2, respectively. In addition, we may use Eq. (4.5) to derive the first derivative of 
A with respect to b 



dA ^f^A ■: 



i=2 



1=1 



(4.11) 



Finally, we need to know the value of Vgff evaluated at the static solution. From the effective 
action (4.6) and for jC^in — 0, we may easily see that Vgff is the effective cosmological 
constant on the brane and it is given by 



7^(4) 3^2 



eff 



4kI 



(4.12) 



Alternatively, the above result may follow by writing = ^eff/A and using Eq. (3.20). 
If we put everything together, the extremization constraint becomes 



A' 



dV,. 



db 



Kg n| 3 A \d(f)J 



0. 



(4.13) 



i=l 



The above expression is indeed zero as can be seen by using the explicit jump conditions 
(3.13)-(3.14). Thus, we may confirm that the solution (3.11), together with (3.5), extremizes 
the radion effective potential. 
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The stability behaviour of those solutions under small, time-dependent perturbations 
will become manifest in the sign of the second derivative of the effective potential. The 
stabilization constraint, that requires a positive second derivative, can be obtained by dif- 
ferentiating Eq. (4.9) one more time with respect to b. If we again use the jump conditions 
(4.10) as well as the scalar field equation in the bulk, we find 



> d'Veff 
db^ 



1=1 

2 



5J 
18 



+ 



1 d'VifdViy 



E 7 y'^^^ 



d''A^ 



\db) ^ dh" 



eff' 



(4.14) 



In the above expression, a term proportional to dVeff/db has been already dropped due 
to the extremization constraint. If we use the jump conditions (3.13)-(3.14), the second 
derivative takes the form 



db'' 



-3 



[y2 sinh(3iJt/2) - Vi smh{3Hyi)f 



+ 2 



/^isinh(3i/|yi|) 
yf cosh(3if|/i) - yj cosh{3Hy2) 



cosh(3i?yi) — cosh(3i?y2) 

l/^cosh(3i?yi) y| cosh(3ify2) 1 



(4.15) 



smh\3Hyi) sinh^(3//y2) 



where 



12//2 



2 4 



d'^V, 



(4.16) 



Kg smh\3Hyi) d(f)'^ 

is the term involving second derivatives of the interaction terms of the scalar field on the 
branes. Let us, for the moment, ignore this term and concentrate on the expression inside 
curly brackets in Eq. (4.15). A simple numerical analysis will reveal the sign of this 
combination. In Figure 1, we display the value of this expression, as a function of H and 
the intcr-brane distance parameterized by the location of the second brane, y2 G {yi, 0) - for 
simplicity, we fix yi = — 1. As one can see, this combination is always negative and becomes 
more negative as both H and the inter- brane distance increases. Eqs. (4. 15)- (4. 16) are valid 
for if^ > 0, that corresponds to a de Sitter four-dimensional spacetime. Solutions with a 
positive cosmological constant on the brane have been derived in the literature before and 
they were shown to be unstable [21]. However, the presence of a bulk scalar field, with 
non-trivial brane interaction terms, may significantly modify this picture: the destabilizing 
behaviour displayed in Fig. 1 may be counterbalanced by the A'-term (4.16) by choosing 
appropriately the interaction terms Vi. If d^Vi/dcfP' are positive and large enough, stable 
brane-world solutions with a positive effective cosmological constant can easily arise in the 
framework of this model. 
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Figure 1: The behaviour of the second derivative of the radion effective potential V^ff is 
displayed as a function of H (for > and the inter-brane distance. 

In the flat-brane limit, if — 0, the expression inside curly brackets reduces to zero. This 
is also in agreement with results in the literature that consider flat-brane solutions which 
are saddle points (flat directions of the potential) [21, 20]. However, even in the flat limit, 
the presence of the scalar fleld allows for a non-trivial X-term and the second derivative is 
written as 



Again, choosing appropriately the interaction terms of the scalar fleld, even the flat-brane 
solutions a la Randall-Sundrum can be made stable. 

Finally, in the case of Anti de Sitter spacetime, = —H^ < 0, the expression for the 
second derivative follows from Eq. (4.15) after making the substitutions: H —>■ iH and 
sinh(3if||/i|) i sin(3if|?/i|). In that case, the expression inside curly brackets turns out to 
be positive-deflnite revealing the existence of a stabilizing force acting on the system. This 
is again in agreement with results presented in the literature [21]. The A'-term in this case 
has the form 

and it may act again as an extra stabilizing force if d'^Vi/dcfP' > 0. It is worth noting that, 
if the above inequality holds, the A'-term acts as a universal stabilizing force independently 
of the value of the effective cosmological constant on the brane. 

4.2 Stability behaviour of the solutions with Vb ^ 

We now turn to the stability analysis of the solutions derived in section 3.2, that are char- 
acterized by a non- vanishing bulk potential Vb- Changing to non-conformal coordinates. 
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the perturbed five-dimensional line-element may be written as 

ds^ = n^{t, y) [ -dt^ + a^{t) dx' dx^ ] + dy^ . 



(4.19) 



The above metric ansatz is similar to that used in Ref. [20], therefore, the details of the 
analysis that lead to the expressions of the stabilization constraints can be found there. 
Here, we present only the main points of our calculation that allow for the generalization 
of those results in the case of a non- vanishing effective cosmological constant. 

For the above ansatz, the scalar curvature may be expressed as in Eq. (4.2) with the 
only difference, relevant to our purposes, being 



1 /12n' 



/2 



6^ \ 



+ 



■in" 
n 



(4.20) 



Substituting the expression of R in the action and performing the same conformal trans- 
formation, the effective potential finally takes the form 



A\hy,,hy2)%f{^1 r^di 

Jbyi 



6 0'2 ^ 



j:nt{V, + K), (4.21) 



i=l 



where now 



A{byi,by2) 



2^4 

^5 



Jbyi 



(4.22) 



The extremization constraint demands the first derivative of the effective potential, with 
respect to 6, to vanish when evaluated at the background solution (3.32). Differentiating 
once with respect to h and using the jump conditions 



Hi 6 



K + Vi{(t)i) 

for i = 1, 2, respectively, we obtain: 



2 dcf) 



(4.23) 



V=Vi 



A' 



dV, 



db 



i=2 



(4.24) 



i=l 



The expression of dA/db is given, in this case, by Eq. (4.11) with nf replaced by nf, while 
Veff is still given by Eq. (4.12). Finally, we may use the relation 



VB(yi)^3E--<l>f^3E 



3 fdA^Y 

'4 J 



(4.25) 



If we put everything together, the extremization constraint becomes 



,dZ 



eff 



db 



i=2 



= 0. 



(4.26) 
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We can easily check that the second equahty indeed holds if we use the jump conditions 
(3.36)-(3.37). 

After differentiating twice Eq. (4.21) with respect to b, we arrive at the stabilization 
constraint. Once again, using the jump conditions (4.23) and the equation of motion of the 
scalar field in the bulk (3.27), together with Eq. (3.26), we arrive at 



i=l 



1 2 fdVA' 



4f=i 



dAy d^A 
'dbj ^ ~d^ 



eff 



(4.27) 



Let us first concentrate on the case with u;^ > 0. By using the corresponding solution 
for the warp factor from Eq. (3.32) and the jump conditions (3.36)- (3.37), we obtain the 
final expression 

d^Z,r .-f/ 12.^ 54c')y+16Uc^--^)z], (4.28) 



db^ Ac| \ \sinh^(a;7/i) 

where we have defined X, y and Z as: 



sinh {ouyi 



i=l 



d^ 
d(f)^ ' 



(4.29) 



y 



sinh (2a; 1^1 1) - |/| sinh(2(x;|j/2|) 
sinh^(a;yi) 



(4.30) 



y2 sinh^(cjy2) - Vi smh^{ujyi 



sinh (cut/i) 2uj{yi - 1/2) + sinh(2cj|yi|) - sinh(2ci;|?/2|) 



(4.31) 



In Eq. (4.28), the constraint relating x ^ that appears in Eq. (3.32), may be 
rewritten as 



4c' 



sinh (coyi) 



(4.32) 



and has been used to eliminate one parameter from the expression of the second derivative. 
The value of this quantity still depends on uj and c^, that parametrize the size of the bulk 
quantities, and Cb{4>), and the inter-brane distance, parametrized by the values of 
For simplicity, we may fix the locations of the two branes at i/j = — 1 and 1/2 = —0.5, while 
the singularity lies at y = 0. As was pointed out in section 3.2, the parameter may take 
both positive and negative values; this may lead to different behaviours of the background 
solution under time-dependent perturbations, as we will now see. In Figure 2, we display 
the value of the second derivative of the effective potential - the expression inside curly 
brackets in Eq. (4.28) - for both signs of and we comment on the results below: 
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Figure 2: The two plots depict the second derivative of the radion effective potential Ve/f 
as a function of u (with u'^ > 0) and for positive and negative values of , respectively. 

• (? > : The second derivative assumes positive values for large enough but small 
uj. In other words, the background solution is more stable when the kinetic term of 
the scalar field is larger and the bulk cosmological constant is smaller. Eq. (4.32) 
then reveals that this particular regime of parameters corresponds to solutions with 
a negative cosmological constant on the brane. 

• < : Positive values for the second derivative can be achieved for large enough 
values of both and u. In this case, the effective cosmological constant on the 
brane is positive by definition and therefore physically interesting, stable solutions 
may emerge in this case. 

It is worth noting that, in the above cases, we obtain positive sign for the second derivative 
of the radion potential even if we assume that d'^Vi/dip'^ = 0, contrary to what happened 
in section 4.1. An extra stabilizing force may arise, for > or < 0, if d^Vi/dc!? > 
or d'^Vi/d(f)'^ < 0, respectively. Let us finally note that in the parameter regimes where 
the solution is clearly stable or unstable, the sole effect of the inter-brane distance is to 
change the magnitude of the radion mass squared: larger inter-brane distances imply larger 
(in absolute value) second derivatives. In the intermediate regime, where the solutions 
struggle between stability and instability, stable solutions may arise when the second brane 
is located closer to the first one and away from the singularity. 

In the special case of u = 0, the expression of the second derivative of the effective 
potential may be obtained from Eq. (4.28) by taking the limit — > 0. Then, we find 

j2t7 2 2„,4 J2\r 9 

= E -^^ + 4-2 {4c' - Qx') {\yi\' - \y2\') . (4.33) 

db^ ~t Vi dcj)^ ^Vi 
The second term has a positive value only if > 6x^ > 0. In this case, the first term 
also contributes positively to the second derivative if d^Vi/dcj)^ > 0. In this case, stable 
solutions with a negative cosmological constant on the brane arise. 
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Figure 3: The plot depicts the second derivative of the radion effective potential Vefs a 
function of the uj and (? parameters. 

Finally, we turn to the case with cu^ < 0. The expression of the second derivative may 
again be obtained from Eq. (4.28) by making the substitutions: lo iuj and sinh(co'?/j) — ^ 
isin(d;?/i). Given the new form of the warp factor, the two branes should now be located 
between the two singularities at ?/ = and ujy = —n. We therefore set yi = —n and we 
allow LJ to vary in the interval (0,1) - y2 is placed in between, with the size of the inter-brane 
distance having no effect on the behaviour of the second derivative. In Figure 3, we display 
the value of the second derivative as a function of and a). We may easily see that stable 
solutions arise for large enough, positive values of c^, with the solutions becoming more and 
more stable as the value of u increases. An unstable regime appears when oj —>■ 1, which 
is equivalent to placing the first brane very close to the singularity a.t ujy = —n. By using 
Eq. (4.32), we may see that the stable regime corresponds to solutions with a negative 
cosmological constant on the brane. 



5 Conclusions 

There has been a lot of work related to the stabilization of the radion field that parametrizes 
the size of the extra dimension in the context of brane-world models. Most work has 
focused on the derivation of solutions with a constant radion field, i.e. a static extra 
dimension, with only a few papers investigating whether these solutions correspond to a 
true minimum of the radion effective potential. In this paper, we presented two new, brane- 
world solutions arising in the presence of a bulk scalar field, and studied their stability under 
time-dependent perturbations of the radion field demonstrating the existence of phenomeno- 
logically interesting, stable solutions with a positive cosmological constant on the brane. 

For our analysis, we used a factorizable ansatz for the line-element along the brane. 
Under the assumption that the total energy of each brane is the sum of a constant brane 
tension and the interaction term of a time- independent, bulk scalar field, the factorization 
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of the 3D scale factor was shown to be directly related to the stabilization of the extra 
dimension. This relation holds both on and off the brane, as one can see by using the scale 
factor jump conditions and the off-diagonal component of Einstein's equations, respectively. 

The first class of solutions presented here corresponds to a vanishing bulk potential for 
the scalar field and a vanishing bulk cosmological constant. Surprisingly enough, brane- 
world solutions with non-trivial warping along the extra dimension did emerge, with the 
'warping' parameter being the expansion rate on the branc. These configurations accept 
a variety of time-dependent solutions for the scale factor on the branc, with fiat (k = 0) 
or curved {k = ±1) 3D spacetime and a zero, positive or negative effective cosmological 
constant. All of the above solutions are characterized by the presence of a bulk curvature 
singularity, which inevitably leads to the introduction of a second brane in the model. One 
can show that the jump conditions, imposed on the warp factor and the scalar field, lead to 
the fixing of the inter-brane distance, in terms of the fundamental parameters of the theory, 
and that the conventional Priedmann equation on the brane is successfully recovered. 

In the second class of solutions that we derived, the presence of a non-trivial bulk 

potential and a bulk cosmological constant was restored. In this case, the warping of the 
5D metric was governed by the bulk cosmological constant shifted by a constant quantity 
given in terms of the kinetic and potential energy of the bulk field. The same variety of 
cosmological solutions on the branc emerge here as well. The sign of the kinetic term of the 
scalar field reveals its nature (normal or tachyonic) and affects the behaviour of the bulk 
potential: a normal kinetic term leads to a potential which is unbounded from below, near 
the bulk singularity, whereas a tachyonic kinetic term leads to an infinitely-high potential 
barrier that may be used to shield the singularity in single-brane configurations. Here, we 
introduced a second brane in order to do so, and we demonstrated that, as in the first case, 
the inter-brane distance is fixed and the form of the Friedmann equation on the brane is 
recovered. 

In the second part of our paper, we investigated the stability of our solutions under small, 
time-dependent perturbations of the radion field. We derived the extremization constraints 
for both types of solutions and demonstrated that, as expected, they correspond to extrema 
of the radion effective potential. The stabilization constraints were also derived. These 
revealed the stability behaviour of the solutions and the type of extrema to which they 
correspond, either minima or maxima. In the first class of solutions, brane configurations 
with positive, zero or negative effective cosmological constant were studied and it was shown 
that, in the absence of a scalar field potential, those solutions come out to be local maxima, 
saddle points or minima of the radion effective potential, respectively, in agreement with the 
literature. However, in our case, the presence of an extra term, involving second derivatives 
of the interaction terms of the scalar field on the branes, acts as a universal stabilizing force, 
independent of the sign of the cosmological constant on the brane, as long as the second 
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derivatives are positive. In this way, solutions with a positive cosmological constant on the 
brane may become stable, saddle-point solutions a la Randall-Sundrum may turn to true 
minima, and AdS-type solutions on the brane may be further stabilized. The results for the 
second class of solutions found in this paper are even more interesting: the aforementioned 
term with the derivatives of the interaction terms may still act as an extra stabilizing force - 
upon appropriate choice of the sign of the second derivatives - nevertheless, stable solutions 
can arise even in the case where this term is zero. The parameter regimes that correspond 
to stable solutions are determined by the values of the bulk cosmological constant and 
the kinetic term of the scalar field. It is worth noting that the sign of the latter quantity 
also defines the sign of the effective cosmological constant of the stable solution: a normal 
(positive) kinetic term gives rise to solutions with negative cosmological constant, while a 
tachyonic (negative) one leads to a positive effective cosmological constant. 

We may, therefore, conclude that the introduction of a bulk scalar field, in a brane-world 
model, may successfully lead to a variety of stable solutions, but more importantly, it may 
lead to cosmologically interesting, stable solutions with a positive effective cosmological 
constant - a type of solutions that has been difficult to derive up to now. 
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